Numerical renormalization group study of random transverse Ising 
models in one and two space dimensions 

Yu-Cheng Lin 1 , Naoki Kawashima 2 , Ferenc Igloi 3 ' 4 , Heiko Rieger 1 ' 5 

1 NIC c/o Forschungszentrum Jiilich, 52425 Jiilich, Germany 
2 Department of Physics, Tokyo Metropolitan University, Tokyo 192-0397, Japan 
3 Res. Inst. f. Solid State Physics & Optics, 1525 Budapest, P.O.Box 49, Hungary 

4 Institute for Theoretical Physics, Szeged University, 6720 Szeged, Hungary 
5 FB 10.1 Theoret. Physik, Universitdt d. Saarlandes, 66041 Saarbriicken, Germany 

(Received ) 

The quantum critical behavior and the Griffiths-McCoy singularities of random quantum 
Ising ferromagnets are studied by applying a numerical implementation of the Ma-Dasgupta- 
Hu renormalization group scheme. We check the procedure for the analytically tractable 
one-dimensional case and apply our code to the quasi-one-dimensional double chain. For 
the latter we obtain identical critical exponents as for the simple chain implying the same 
universality class. Then we apply the method to the two-dimensional case for which we get 
estimates for the exponents that are compatible with a recent study in the same spirit. 

§1. Introduction 

The effect of quenched randomness on disordered quantum magnets close to a 
quantum phase transition is much stronger than on classical systems at temperature 
driven phase transitions. As first observed by McCoy 1 ) in a somewhat disguised 
version of a random transverse Ising chain, non-conventional scaling and off-critical 
singularities that lead to divergent susceptibilities even away from the critical point 
now appear to be a generic scenario in any dimension, at least in disordered quantum 
magnets with an Ising symmetry. The reason for this, as pointed out by Fisher 2 ) 
only recently, is a novel fixed point behavior of these systems under renormaliza- 
tion, namely one which is totally determined by the randomness and its geometric 
properties: the so called infinite randomness fixed point 2 }' 3 \ 

Within this scenario the quantum critical behavior of disordered transverse Ising 
models is essentially determined by strongly coupled clusters and their geometric 
properties 2 )' 3 ^. Let L be the linear size of such a cluster. Then it contributes to the 
low energy spectrum with an exponentially small excitation gap of size In AE ~ , 
defining the exponent ip. Moreover, at the critical point, it has a total magnetization 
of size /i ~ defining the exponent (j>. Finally the linear length scale of strongly 
coupled clusters occurring at a distance 5 away from the critical point is £ ~ \5\~' / 
giving rise to a third scaling exponent v. All bulk exponents can now be expressed 
via tp, (ft and u, c.f. flb/v = Xb = d — <fnp, ft yp = v[l — ip) and in the Griffiths phase 
z' oc 5~ u ^ . For the Id case, as treated above, it is ift = 1/2, <j) = (v5 + l)/2 and 
v = 2 for uncorrelated disorder. 

The basic geometric objects, the strongly coupled clusters, still have to be de- 
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fined and this will be done within a renormalization group scheme. However, for 
site or bond dilution it is immediately obvious what these clusters are: simply the 
connected clusters. Hence the critical exponents defined above are directly related 
to the classical percolation exponents 4 ); Let 8 = p — p c be the distance from the 
percolation threshold, f pcrc the exponent for the typical cluster size, -D per c the fractal 
dimension of the percolating cluster, /3 pe rc the exponent for the probability to belong 
to the percolating cluster. Then one has for the critical exponents defined above 

^ = ^pcrc > ifi = -Cperc , <j> = (d /3pcrc/ J/ perc)/-Cperc (1) 

Next we consider the question, what happens for generic disorder (i.e. no dilution, but 
random bonds and/or fields) and we consider the model defined by the Hamiltonian 

H = ~ E J ij°!°j ~ E Mf • (2) 

Here the {erf} are Pauli spin matrices, and the nearest neighbor interactions and 
transverse fields hi are both independent random variables distributed uniformly: 

1, for < Jij < 1 
0, otherwise 
Kq 1 , for < hi < ho 
0, otherwise 

For this case the distance 5 from the critical point is conveniently given by S = 
| In ho. In one space dimension this model has been investigated intensively over the 
recent years 6 ). 5 ). 7 ) " 9 ) ; an d many analytical as well as numerical tools are at hand to 
analyze it. Beyond the simple one-dimensional geometry one has to rely on numerical 
techniques like quantum Monte-Carlo simulations (as in the two-dimensional case 10 )) 
or the numerical implementation of the renormalization group scheme, which we 
outline in the next section. 



p{hi) 



§2. The renormalization-group scheme 

The strategy of the renormalization-group a la Ma, Dasgupta and Hu 11 -* is to 
decrease the number of degrees of freedom and reduce the energy scale by perform- 
ing successive decimation transformation in which the largest element of the set of 
random variables {h^ J^} at each energy scale is eliminated and weaker effective 
couplings are generated by perturbation theory. 

The renormalization-group procedure is as follows: Find the strongest coupling 

I? = maxjjjj, hi} 

in the system. If i? = Jy, then the neighboring transverse fields hi and hj can be 
treated as a perturbation to the term —J^afaj in the Hamiltonian (2); The two 
spins involved are joined together into a spin cluster with an effective transverse 
field 

~ hihj 

n (ij) ~ ~r~ 
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Fig. 1. Schematic of renormalization-group decimation for spin chains 

and an effective magnetic moment 

H%3) = + H ■ 
The bonds of the new cluster with other clusters o~k are 

J(ij)k ~ max(Jj fe , Jj k ) . 

If instead Q = hj, then the associated spin o~j is eliminated and effective bonds 
between each pair of its neighboring spins are generated by second-order perturbation 
theory. The strength of the effective bonds for each pair (i, k) is 



J ik w max( Jj 



ik ■ 



hj 



where the Jik are the bonds that may have already been present. This procedure is 
sketched for the Id case in Fig. 1. We continue the procedure until there is only one 
remaining spin cluster. 

At each stage of the RG, an effective field (bond) is a ratio of a product of 
some number / of original fields (bonds)to a product of original / — 1 bonds (fields). 
The / grows under renormalization at criticality. As a result, the log- field and 
log-bond distributions R^(lnh) and Pft(ln J) become broader and broader under 
renormalization as the critical point is approached. This increasing width of the field 
and bond distributions reduces the errors made by the second-order perturbation 
approximation. The RG becomes thereby asymptotically exact. 



§3. The one-dimensional case 

The RG can be carried out analytically in one space dimension 6 ), therefore 
we can use the Id case with periodic boundary conditions as a simple check for 
our numerical implementation. In Fig. 3 we show the probability distribution of the 
logarithm of the last remaining cluster field at the critical point ho = 1, which scales, 
according to Fig. 2, like L -1 / 2 , where L is the system size. From this one concludes 
that the exponent tp defined in the introduction, is given by tp = 1/2. Inspecting the 
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Fig. 2. Left: Distribution of the logarithmic strength of the last remaining cluster fields, r = 
ln(-||^) (J?o denotes the energy scale of the original Hamiltonian) . The distribution gets broader 
on a logarithmic scale with increasing system size, indicating a infinite dynamical exponent z. 
The data is obtained from 100 000 samples for each system size. Right: Scaling of the data in 
the left figure, assuming the exponential scaling form obtained from the analytical work. 6 - ) 
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Fig. 3. Scaling of the number of active spins 
(proportional to average magnetic moment 
per spin, n) in the last remaining spin clus- 
ter at the critical point. We find fj, ~ L 81 
implying (j> w 1.62. 
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Fig. 4. In the ordered phase (5 < 0), the spon- 
taneous magnetization scales as Mo = |<5| fl 
with (3 w 0.39. This numerical estimate of 
f3 is in agreement with the analytical pre- 
diction: ft — 3 ~ 2 V ^ . Our data is obtained 
by averaging over 100 000 samples of size 
L = 1024. 



number of active spins in the last remaining cluster at the critical point we obtain 
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Fig. 5. Left: Distribution of logarithmic effective fields fih of the last spin cluster in the disordered 
phase at ho = 2. The curves for different sizes look very similar but shifted horizontally related 
to each other. Right: Scaling plot of the data in the left figure. The typical fields and spacing 
of rare large strongly coupled clusters in the disordered phase is related via Q ~ L z . The 
dynamical exponent z'(h) is obtained from the asymptotic form lnP(ln .!?;,) = l/z'(h)(ln Qh) + 
const. 5) 



1/z' 




Fig. 6. The value of 1/z' (5) for d = 1 obtained 
from the probability distribution of loga- 
rithmic effective fields of the last remain- 
ing clusters in the Griffiths-McCoy region. 
The solid line, which fits well, is a plot of 
the exact 2'(5)-relation reported in Ref. 8 ^: 
2'log(l - z'~ 2 ) = -2h . 



the size dependence fj, ~ L ' 81 from Fig. 3, and thus 4> ~ 1.62. 

In the Griffiths phase ho ^ 1, the probability distribution of the energy gap f2 
still has an algebraic singularity at Q = 0, and its finite size scaling behavior is 

QP L {Q) oc L d Q d / z '^ = (L z '^f2) d / Z '^ (3) 

where d is the space dimension (in this section d = 1) and z'(5) a generalized dynam- 
ical exponent that varies continuously with the distance 5 from the critical point. 
This exponent parameterizes the strength of all singularities in the off-critical region 
5^0, for instance in the disordered phase 5 > one has for the imaginary time 
autocorrelations Gi oc (r) = [(of (T)o"j(0))T=o]av ~ r -1 / 2 '^, for the local susceptibil- 
ity xioc ~ T l l z ( <5 )~ 1 , for the specific heat C ~ T l / Z ^ and for the magnetization 
in a longitudinal field M ~ H l l z The most convenient way to determine this 
exponent is, however, via the distribution Pl{Q). At the critical point this distri- 
bution has to merge with the critical distribution discussed above — and therefore 
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lim^o z '{$) = 00 • Using this finite-size scaling form for the distribution of the last 
bonds/fields in the RG procedure we can extract the dynamical exponent as is done 
in Fig. 5. 

In the ordered phase ho < 1 the distribution of fields and bonds are related to 
the distribution in the disordered phase ho > 1 via duality, see Fig. 7. 




ft, ln(L 2 ' ft,) 



Fig. 7. Distribution of last log-bonds (left) in the ordered phase at ho = 0.5 for d = 1 and its 
scaling plot (right). One observes that the scaling behavior of bonds in the ordered phase and 
that of fields in the disordered phase (see Fig. 5 ) are related through duality. 



§4. The double chain 

The RG scheme for double chains with some new elements (compared to the Id 
case treated above) is depicted in Fig. 8. 
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Fig. 8. Schematic of renormalization-group decimation for double chains used in the numerical 
simulations. 

As in Id, we observe that the log-field and log-bond distributions get broader 
with increasing system size at criticality. To estimate the critical point we compute 
the field distribution at the last stage of the RG varying the initial transverse field 
ho- We estimate the critical point to be at ho = 1.9, beyond which the broaden- 
ing of the log-field distribution appears to be saturating, as for Id in the Griffiths 
phase. Moreover at h c = 1.9 the log-field and the log-bond distributions become 
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asymptotically identical except for a constant multiplicative factor that reflects the 
short-ranged non- universal physics (See Fig. 9). This is obvious in the single chain, 
where it follows from the self-duality of the simple chain at the critical point. How- 
ever, the double chain is not self-dual, nevertheless the scaling forms of the two 
distributions become identical at the critical point. We speculate that this remains 
true also in the two-dimensional case to be discussed below. 
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Fig. 9. The log-field (left) and log-bond (right) distribution for the double-chain at h c = 1.9. The 
data is obtained from 25 600 samples for each system size. 
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Fig. 10. Scaling plots of the distribution of log-fields (left) and of log-bonds (right) for the double- 
chain at the critical point (h c = 1.9). The data scales quite well with the same form -Th(j) = 

W «2_) ~ VI as in Id. 
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The scaling of the critical distributions depicted in Fig. 9 yields the critical 
exponent ip = 0.5, as shown in Fig. 10. This is the same as for the simple chain. In 
addition, for the average magnetic moment of the last remaining cluster at h c = 1.9, 
we find the same system size dependence for the double chain as for the Id case, i.e. 
the same critical exponent (j), see Fig. 11. This implies that the double chain and 
the simple chain belong to the same universality class. 

In the Griffiths phase ho > h c = 1.9 we extracted the generalized exponent 
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Fig. 11. A log- log scaling plot of the average 
magnetic moment per cluster, /i, with the 
linear system size L at the estimated criti- 
cal point h c — 1.9 for the double chain. We 
find the same exponent <j) as in lei. 



0.8 



0.6 



■S 0.4 
0.2 



1.5 



Fig. 12. Estimates for 1/z' against ho for the 
double chain. 
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z'(ho), which is depicted in Fig. 12. Close to the critical point h c we observe the 
same linear dependence of 1/z' (8) on the distance 8 = ho — h c from the critical point 
as in Id. Since for 8 <C 1 one expects z'(ho) oc 8~^ v this implies that v = 2 the same 
as the simple chain. 

§5. The square lattice (2d) 



Next we present our preliminary results for the two-dimensional (2d) case with 
periodic boundary conditions, where we, in contrast to the treatment in Motrunich 
et al. 3 ), keep all bonds generated during renormalization. The RG scheme for the 
2d case is very similar to the one for the double chain and is depicted in Fig. 13. 

In comparison to Id and the double chain models, the location of the critical 
point cannot be fixed precisely for two dimensions according to our numerical obser- 
vation so far. We obtain a critical field approximately at ho = 5.3 by applying the 
criterion that field and bond distribution should have the similar scaling form (as for 
the Id case and the double chain) . The scaling of the last log-field distribution yields 
ip « 0.5 and the scaling plot of the number of the active spins in the last remaining 
cluster yields <j> ss 2.0 and fi ~ L im . 



Numerical renormalization group study of random transverse Ising models 9 









' i 1 






h> 


} 






J,* 


hi 



n = Jh 



n = ^ 





r hi hi 

Fig. 13. Renormalization-group decimations for the two-dimensional (square) lattice used in the 
numerical simulations. 
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Fig. 14. Scaling of the log-fields for d — 2 at 
the last stage of the RG at what we esti- 
mate to be the critical point. 
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Fig. 15. A scaling plot of the number of the 
active spins in the last remaining cluster at 
our candidate critical point ho = 5.3. 
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Our preliminary results for the two-dimensional case agree with those obtained 
recently by Motrunich et. al. 3 ) and with those obtained by us via quantum Monte- 
Carlo simulation 10 ) . 
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